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Abstract 

We present 0{a s ) results on the decays of polarized and Z bosons into massive 
quark pairs. The NLO QCD corrections to the polarized decay functions are given 
up to the second order in the quark mass expansion. We find a surprisingly strong 
dependence of the NLO polarized decay functions on finite quark mass effects even 
at the relatively large mass scale of the and Z bosons. As a main application 
we consider the decay t — y b -\- W~^~ involving the helicity fractions Pmm 

of the W + 

boson followed by the polarized decay Vk + (t) — > <?i<?2 for which we determine the 
0{a s ) polar angle decay distribution. We also discuss NLO polarization effects in 
the production/decay process e + e~ —tZ(f)—} qq. 



1 Introduction 



The polarization of W and Z bosons produced in electroweak production processes is in 
general highly nontrivial. One therefore has a rich phenomenology of polarization effects 
in (W, Z) production and decay which will be explored in present and future experiments. 
The polarization of the W and Z bosons can be probed by decay correlations involving 
the decay products of the polarized (W, Z) bosons. A widely discussed prominent example 
of such decay correlations is the decay t — > b + W + followed by W + — > i^v^ where the 
decay W + — » is used to analyze the helicity fractions of the W + resulting from the 
decay t — > b + W + (see e.g. Refs. [HE])- It would be interesting to explore the possibility 
to also make use of the quark-antiquark decay modes W ± — > q\qi and Z — > qq to analyze 
the polarization of the (W , Z) bosons, in particular using the tagging modes W + — > cb, cs 
and Z — y cc, bb involving heavy quarks. 

This paper is devoted to the calculation of NLO QCD effects in the decay of polarized 
(W, Z) gauge bosons into massless and massive quark-antiquark pairs. In order to provide 
a quick access to the importance of quark mass effects we have made use of a quark 
mass expansion of the rather lengthy fully analytic NLO results listed in [3j. We thereby 
demonstrate that, in polarized gauge boson decays involving charm and bottom quarks, 
the NLO finite mass effects are non-negligible. The reason is that the NLO finite mass 
corrections in polarized decays set in at linear order with rather large coefficients, contrary 
to the case of unpolarized decay where the finite mass effects set in only at 0{m 2 /m^ vz ). 
Depending on the particular polarized decay function, the mass corrections can become as 
large as the leading term of the NLO mass expansion for decays involving b and c quarks. 

NLO and finite mass effects will affect the decay correlations between the momenta of 
the production and decay process. As specific examples of how such correlations are affected 
by NLO and finite mass effects we consider the cascade decay process t — > b J r W + { J \){— > cb) 
and the production/decay process e + e~ — > Z{\) — > bb,cc. 
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2 Angular decay distribution 

The polar angle decay distribution of a polarized (W, Z) boson decaying into a fermion- 
antifermion pair (quark or lepton pair) is given by 

W(9) = Pmmd 1 mm ,{6)d l mml {6) H m i m i 

m,m'=0,±l 

= |(1 + cos 2 6) ((p ++ + p__) (H ++ + #__) + 2p 00 ^oo) 
+^ C os6(( P++ -p^)(H ++ -#__)) 

+ - sin 2 9 ((p++ + P—)H 00 + p m (H ++ + H '__ - floo)) • (1) 
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The diagonal spin density matrix elements of the polarized (W, Z) boson are denoted 
by PmmJj They are defined in a (x,y, z) coordinate system associated with the production 
process while the polarized decay structure functions (for short: polarized decay functions) 
H m i m i are defined in a (x',y, z') coordinate system associated with the decay process. The 
two coordinate systems are rotated into each other by a rotation around the y axis by the 
polar angle 9. For example, for the sequential decay t — > b + W + ( J \)(— > cb) the z axis could 
be chosen to lie along the momentum of the W + in the top quark rest frame and the z' 
axis could be chosen to lie along the quark direction in the W + rest frame. 

The polarized decay functions H±± and ifoo stand for the probability of the decay 
of a polarized vector boson (W, Z)(m') into a fermion-antifermion pair where the vector 
boson has the spin quantum numbers m' = 0, ± in the rotated {x' , y, z') coordinate system. 
The spins of the fermion-antifermion pair are summed over when calculating the decay 
probability. 

We take the spin density matrix to be normalized, i.e. p ++ + p o + P = 1- It 

is also convenient to define a normalized angular decay distribution W{9) = W(6)/W 
where W = f\ W(6)dcos9 = J2m H mm such that f\ W(6)dcos9 = 1 (later on we denote 



the literature the diagonal elements of the normalized spin density matrix of the (W, Z) boson p n 
are frequently referred to as the helicity fractions T m of the gauge boson. 



W = Em Hmm by H u+L ). Correspondingly we define normalized polarized decay functions 

Hm'm' H m i m i / E m H mm . 

The distribution Eq. (CQ) or its normalized form W{9) is a second order equation in 
cos#, i.e. it is the equation of a parabola with coefficients given by sums of the products 
p mm H m i m i. The parabola is upward bent for (1 — 3p o)(l — 3i/ o) > and downward bent 
for (1 — 3poo)(l ~ 3i7oo) < 0. As a measure of the flatness of the decay distribution we 
define a convexity parameter c/ given by the differential change of slope (or the second 
derivative) of the decay distribution. From Eq. (fl]) one has 

As a second global measure we introduce the forward-backward asymmetry of the decay 
distribution defined by 

W(F) — W(B) 3, ... ~ , , Q , 

= w\f ) + w\b) = 4 ( ^+ - ~ (3) 

where W(F) = W(0 <9< tt/2) and W(B) = W(n/2 <6<tt). 

Of interest is also the location of the extremum of the parabola in Eq. (fl|). The ex- 
tremum is located at 



cos 9 



(P++-P-) (#++-#-) (4) 



extr (1 _ 3poo ) (1 _ 3Hm 



The three measures are not independent for the normalized parabolic decay distributions 
W(9). They are related by 

cos 9 



AfB (5) 



Note that all three measures factor into a production part described by the density matrix 
elements p mm and a decay part given in terms of the polarized decay functions H mm . 

As a well-known illustration consider the cascade decay £ — > b + W + ('\) followed by 
VT + (t) — > /i/2- At the Born term level and for massless fermions (quarks or leptons) 
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the only nonvanishing polarized decay function in the Standard Model (SM) is H ++ if one 
chooses z' to lie along f 2 G {£ + , (72 }• in the normalized form one has H++ = 1. One obtains 

W(6) = I (1 + cos Of p ++ + ~ (1 - cos Of p__ + ~ sin 2 9 p 00 . (6) 
o o 4 

The decay distribution fl6]) corresponds to a normalized parabola with convexity parameter 
Cf = 3(1 — 3/?oo)/4, a maximum at cos^ | extr = — (p++ — p )/(l — 3poo) an d a forward- 
backward asymmetry of A FB = |(p++ — p ). 

It is evident that an analysis of an angular decay distribution such as in Eq. (JSJ) can 
be used to experimentally extract information on the helicity fractions p mm (m = ±1,0) 
of the W boson. Such an analysis has been widely applied using the leptonic decay modes 
of the W boson iy + (|) ->■ P + ve (see e.g. Refs. [HE])- Clearly the decay distribution (J6]) 
will be affected by radiative corrections and finite mass effects in as much as the three 
normalized decay structure functions H mm will change from the simple pattern H ++ = 1; 

H Q0 = H = on which Eq. ([6]) is based. The modification of Eq. ([6]) in the quark sector 

through radiative corrections and finite mass effects is the subject of this paper. In this 
context it is important to note that for t — > b + W + ('\)(— > qiq?) the radiative corrections 
to t -> b + W + (-\) and W + (\) ->■ q x q 2 factorize at NLO in QCD (but not in higher orders), 
i.e. there is no NLO cross-talk between the production and decay processes [I]. 

Note that for unpolarized (W, Z) decay, when p ++ = p 00 = p = 1/3, Eq. (CD) leads to 

a flat decay distribution 

W(B) = ±(H ++ + H Q0 + H„). (7) 

In the next two sections we shall first write out the diagonal spin density matrix elements 
Pmm of the (W, Z) boson in two prominent sample production processes and then, at a later 
stage, proceed to calculate the relevant set of decay structure functions H mm at 0(a s ). The 
results are then combined to present analytical and numerical results for the respective 
0(a s ) decay distributions. 
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3 Spin density matrix of the W boson 
in the decay t ^ 6 + Vi/ + (t) 

Let us start by discussing the spin density matrix of the W + in the decay t — > b + W + (~f). 
As indicated by the notation, the W boson emerges in a polarized state in this decay. 

The spin density matrix elements of the W + in t — > b + W + have been well studied. 
We take the z axis to lie along the momentum of the W + in the top quark rest frame. At 
LO one has [7] 

p ++ (Born) = -»■ 0.0007, (8) 

Poo(Born) = 1 —— = 0.696 0.6887, 

1 + 2x 2 

2x 2 

p(Born) = r = 0.304 0.3106, 

1 + 2x l 

where x = mw/m t . For the numerical values we use the central values of my/ = 80.399 ± 

0.025 GeV and m t = 172.0 ± 0.9 ± 1.3 GeV provided by the Particle Data Group ©. In 

Eq. © we have also given the NLO QCD results indicated by arrows (cf. Refs. [H [101 

[TT])o The radiative correction to p ++ can be seen to be very small. The absolute (relative) 

corrections to p o an d p amount to —0.73% (—1.05%) and 0.66% (2.17%). 

The spin density matrix elements in Eq. (JSJ) are calculated for unpolarized top decays. 

If the decaying top quark is polarized, the spin density matrix elements will depend on the 

orientation 9p and the degree of polarization P t = \P t \ of the top quark. The dependence 

is very simple for the m.& = Born term case where one has 

P „ P 1 + Pt cos dp p I - Pt cos dp 

P + + = 0, Poo = Poo D[dp) , P„ = P- D(6 P ) ' (9) 

and where the denominator in Eq. ([9]) is given by 

D(9 P ) = Poo (l + P t cos9 P ) + p_4l - P t cos9 P ). (10) 



2 The NNLO corrections to the spin density matrix elements of the W + have recently been calculated 
in Ref. [U. 



It is clear that the relative weight of the two helicity fractions poo and p can be changed 

by appropiately tuning the polarization of the top quark. 

In the general case both helicity states of the polarized top quark are involved. This 
case is somewhat more complicated and has been worked out in Refs. jH [9] . 



4 Spin density matrix of the Z boson 
in the production process e + e~ — > Z^[) 

Next we discuss the polarization of the Z boson produced in e + e~ annihilation where the 
polarization density matrix of the Z boson is purely transverse in the e + e~ system. We 
take the z axis to lie along the e~ beam (z || e~). The Born term matrix element is given 
by (see e.g. Ref. [5]) 

Mim) = -i 9 -fv{e + ) {v elll - a el ^ 5 ) u( e -)e*"(m) (11) 

(m = 0, ±) where, in the SM, one has 

vi = — 1 + 4 sin 2 ®w ai = — 1 for £—e,fj,,T, (12) 

and where Qw is the Weinberg angle. In our numerical analysis we take sin 2 Qw — 0.231 [6]. 
One can then work out the normalized diagonal spin density matrix elements of the Z boson 
in e + e~ — > Z(f) which are given by 

(v e — a e ) 2 (v P + dp) 2 . , 

p++ = A/^ = °- 4244 ' poo = o> p- = knrr^k = °- 5756 - 13 

2{v 2 e + al) 2(t» e 2 + of) 

Note that the two transverse density matrix elements in Eq. (flBl) are approximately equal 
due to the smallness of the vector current coupling constant v e = —0.076. 

Similar to the case t — > b + W + , the spin density matrix elements of the Z boson can 
be tuned by polarizing the e + e~ beams. Take, for example, longitudinally polarized beams 
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and denote the longitudinal polarization of the e T beams by h T , where the polarization is 
measured w.r.t. the momenta of the e T beams. One then has 

P (i + /Q(i-fr+) p p _ p (i-h-)(i + h+) 

where 

D(h~, h + ) = p++(l + hr){l - h+) + - /i~)(l + (15) 

For example, with a 100% longitudinally polarized electron beam one obtains p p __ = 1 for 
h~ = — 1 and = 1 for h~ = +1. 



5 Polarized W decays into massive quark pairs 

We first treat the case W + —> q±q2- The LO Born term amplitude is given by (see e.g. 0) 

Mw {m) = ~i%V qiq2 uM^^^vM ef(m), (16) 

where gw = e/ sinQw is the electroweak coupling constant and the V qiq2 are Kobayashi- 
Maskawa matrix elements. Let us define a reduced matrix element Ai(m) by splitting off 
the coupling factors and the factor 1/2 from the chiral projector such that 

M w (m) = uM^il - l5 )v 2 ( V2 ) e^(m). (17) 

The LO polarized decay functions H mm are then obtained from 

Hmm = N c Y, M w {m)M ] w {m). (18) 

quark spins 

Our aim is to calculate the NLO polarized decay functions H mm of a W + boson decaying 
into a heavy quark pair where the W + has definite spin quantum numbers m — ±, 0. We 
first discuss a coordinate system where the z axis lies along the quark momentum (z || q\) 
(system I). The Born term and the a s contributions to the polarized decay functions H±± 
and ifoo and the unpolarized total decay function Hu + l = H ++ + H 00 + H are expanded 
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up to the second order in the quark mass ratios y/jHp = rni^/rriw where we make use of 
the unexpanded analytical results given in Ref. [3]. One obtains 



H 



++ 



H, 



oo 



8iV c g 2 



8iV c g 2 



+ . 



+^ (l + (tt 2 - 16) JJT X + (5 + 2tt 2 /3) fii + fi 2 ~ 2/xi In Mi + 



67T 

o+ (/il + /i2) +... 



(19) 



+ ^U- 2tt 2 ^ - (39 + 8vr 2 /3) m + /x 2 - 30/xi In 



+6/i 2 In /U2 - 2/xi In //! + ... 



(20) 



iJ 1 



8iV c g 2 



1 - //! - fJL 2 + ■ ■ ■ 



o:. s 



+^( 1 + (tt 2 + 16)^ + (49 + 2vr 2 ) /ix + 13//2 + 14/xi In ^ 
-24yU 2 In /i 2 + 2/xi In 2 /ii + . . . 



(21) 

We truncate the mass expansion at 0(/ij) since third order quark mass effects can be 
expected to be quite small judging from the fact that /i 3 / 2 = (1.5/80.399) 3 = 6.49 ■ 10 -6 
and /4 /2 = (4.8/80.399) 3 = 0.21 • 10~ 3 . 

For the sum of the three polarized decay functions J2 m H m m '■= Hjj + l one has 



Hi 



U+L 



8N c q 



U+L 
2 



H l ++ + Hl + H l _ 
{lii + fi 2 ) 



+ • 



+ -^-(6 + 15/ii + 15/i 2 - 18/iiln/xi - 18/i 2 ln/i 2 + . . . ) 

67T v 



(22) 



Note that the sum of the polarized decay functions is independent of the choice of the z 
axis as is indicated in Eq. ( 1221 . 

The mass corrections set in quadratically in the LO polarized decay functions and also 
in the radiatively corrected unpolarized decay function H u+L . In contrast to this, the 
mass corrections to the radiatively corrected polarized decay functions H±± and H Q0 set 
in linearly. Surprisingly, some of the linear mass corrections carry rather large coefficients 
such as the coefficient (tt 2 + 16) = 25.87 multiplying the linear mass term JJx[ in H l __ in 



Eq. (1211) . Contrary to naive expectations one therefore needs to keep finite mass effects 
in the radiatively corrected polarized decay functions for massive quark pair production 
even at the relatively large mass scale of the W mass. Note that the radiatively corrected 
unpolarized decay function Hu + l is symmetric in the quark masses whereas the polarized 
decay functions show a large quark mass asymmetry at NLO. 

In order to make contact with the unpolarized decay rate r(ll /+ — > qiq~2), we define 
polarized decay rates 

Trnrn — Q^ m 2 9w 1^9192 1 H mm (23) 

(|p| = *iF(l + \A + ^2 ~ — 2/i 2 — 2/i 1 /i 2 ) 1 / 2 ) which, for the unpolarized rate, gives 



T (w + -> qi q 2 ) = i(r ++ + r 00 + r__) = ^M-^ g 2 w \v qiq2 \ 2 h u+l . (24) 



The quark mass corrections will be most important for the decay W + — > cb. For the 
quark masses we take m c = 1.5 GeV and = 4.8 GeV. For the W + — > cb unpolarized 
decay function H u+L one obtains 

H u+L = 8N c q 2 \ 0.998 + . . . + ^ (6 + 0.470 + ...)!. (25) 

In order to highlight the numerical importance of quark mass effects we have separately 
listed the leading and the 0(/i,) quark mass effects in the NLO terms in Eq. (f2"5]) . Even 
though the mass corrections to the sum of the polarized decay functions set in only quadrat- 
ically, the mass effects in the NLO corrections can be seen to amount to a non-negligible 
0(8%) where the largest contribution comes from the bottom quark term — 18/i2hi/i2 = 
0.0071. Using a s (m^) = 0.117 one finds an overall increase of the zero mass Born term 
decay rate by 3.8%, i.e. one has 

r( NLO; 0{ni)) = 1.038 • T( Born; ^ = 0) , (26) 

where the bulk of the increase comes from the NLO zero mass term. 
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For the W + — > cb polarized decay functions one obtains 



H l ++ = 8iV c g 2 [0+... + ^(1-0.101 + ...)], 
H l m = 8N c q 2 [0.002+ ... + ^(4 -0.469 + ...)], 
H l __ = 8N c q 2 [0.996+ ... + -^(1 + 1.040 + ...)]. (27) 

The mass corrections to the NLO polarized decay functions can be seen to be large. The 
largest mass correction occurs for H l __ which is the only polarized decay function which 
is nonzero at the m q = Born term level. The mass correction to the leading NLO 
contribution in H l __ is of O(100%). 

The large NLO mass correction to H l __ does not, however, feed through to the nor- 
malized decay distribution which is governed by the normalized polarized decay functions 
B. I mm = H^ nm /H u+L . This can be appreciated by writing out the polarized decay functions 
in a generic notation where the (small) mass corrections to the Born term contributions 
are neglected. One has 

H ++ = 87V c? 2 (0 + ^(l + /i ++ )) 
#oo = 8iy c g 2 (0 + ^(4 + /i 00 )) 

= 8iV c g 2 (l + -^(l + / i„_)) (28) 

D7T 

where the fi m m denote the respective NLO finite mass corrections. One then expands the 
normalized polarized decay functions H m m in the strong coupling constant a s . In this 
approximation one has 





= o + ^(i + ^ ++ ) 






= + ^(4 + ^) 

07T 






= 1 + ^(- 5 -^++ -/^oo)- 


(29) 



It is apparent that the NLO finite mass corrections to the normalized polarized decay 
functions H mm are solely determined by the O(10%) mass corrections //++ and /ioo- Since 

11 



//++ and /ioo are negative, the NLO mass corrections are destructive. Numerically the NLO 
mass corrections amount to only 0{— 10%) of the leading NLO mass term (see Eq. (T27D ). 

The 0(100%) NLO mass correction to H drops out when normalizing the polarized 

decay functions. 

We are now in the position to write down numerical results for the angular decay 
distribution ([[]). We take into account 0(a s ) results both for the density matrix elements 
p mm and the decay structure functions H mm . At 0(a s ) one thus has to take the sum 
Pmm(Born)H m ' m i(a s ) + p mm (a s )H m i m >(Born). As mentioned before, there is no 0(a s ) 
cross-talk between the production and the decay process because the intermediate gauge 
boson is colour neutral As before we concentrate on the decay W + — > cb. For system I 
one finds (ct s (m^/) = 0.117) 

_ q f 0.305 1 q f 0.001 1 q f 0.694 i 

W l (6) = -(l + cos6) 2 \ 0.318 +-(l-cos#) 2 0.019 +-sin 2 # 0.663 (30) 
8 I 0.318 > I 0.018 > I 0.664 > 

For the sake of comparison we have listed three numerical values each for the angular 

coefficients. The top entry is for (Born; /ij ^ 0), the middle entry is for (0(a s ); /ij = 0), 

and the bottom entry is for (0(a s ); /ij 7^ 0). The same three-tiered notation will be used 

in subsequent formulas and in Table [TJ 

Eq. f )30|) describes a downward bent parabola with unit area and, as Fig. [T] shows, a 

maximum slightly displaced to the right of cos# = 0. By comparing the respective numbers 

in the normalized decay distribution W(9) in Eq. ( )30|) . NLO quark mass effects can be seen 

to be almost negligible even if they are important for the polarized decay functions. More 

important are the radiative corrections which lead to a 2.0% enhancement (on an absolute 

scale) of the normalized decay distribution at the forward point (cos8 = 1) and a 2.7% 

enhancement at the backward point (cos8 = —1). Midways at cos# = one finds a 1.2% 

depletion of the decay distribution. This is illustrated in Fig. [1] where we compare the 0(/i^) 

normalized angular decay distributions for the Born term case and the 0(a s ) case. The 

net effect of the radiative corrections is to make the normalized angular decay distribution 
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COS 

Figure 1: Normalized angular decay distribution W{6) = W(0)/W for W + — > cb at LO (full 
line) and NLO (dotted line) in system I including 0{m 2 c ^) finite quark mass contributions. 

flatter with little dependence on quark mass effects. 

In Table [1] we present our numerical results for the three global parameters c/, Afb and 



cos 9 



that characterize the decay distribution. In Table 1 we use the same three-tiered 



notation as in Eq. (1301) . The numerical results for system I (z \\ c) are listed in column 1. 
We emphasize that the inverse of the denominators factors occurring in the calculation 
of the global measures in Eqs. (J2H1J have been left unexpanded in a s when calculating 
the entries in Table 1. The initial and final state QCD corrections can be seen to reduce 
the convexity parameter c/ and the forward-backward asymmetry by 8.37% and 1.25%, 
respectively. The maximum of the decay distribution is shifted to the right by 7.7%. By 
comparing the numbers in tier 3 with those in tier 2 one can see that the bulk of these 
shifts come from the NLO zero mass corrections. NLO finite mass effects are small and 
tend to slightly reduce the NLO zero mass corrections. 

The largest NLO corrections come from the final state corrections. One can therefore 
obtain a rough understanding of the numbers in Table 1 by neglecting the initial state 
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corrections and the (small) finite mass effects in the Born term contributions. Expanding 
the relevant contributions to 0(a s ) one obtains 

4 = 2(1- 3/^(1-^(12 -1.41)) = "0.762, (31) 

o 

A FB = —(P++-P-) (1-^(6- 0.671)) = 0.221, (32) 
cos # 1 = P++ ~ p — (l + ^i(6 - 0.736)) = 0.289 . (33) 

max 1 — 3p 0Q V 67T 7 

Looking at the different contributions in Eq. (T3TM331) one can understand the main features 

of the numerical results listed in Table 1 that were already discussed above. 

In Eqs. f l3TH33|) we also quote approximate numbers which are calculated from the finite 

mass corrections p ++ and /i listed in Eq. (j2"7|) and the Born term values for the spin 

density elements p mm listed in Eq. (JSJ). The approximate numbers listed in Eqs. (|3lT[33|) 

can be seen to deviate by small amounts from the exact numbers listed in Table 1 where 

i 

the largest deviation occurs for A FB and cos 6 

max 

Summarizing our results in system I one finds that NLO and NLO quark mass effects 
are important for the polarized decay functions. The NLO quark mass effects, however, do 
not feed through to the angular decay distribution which is mainly affected by the leading 
mass term in the final state radiative correction. One concludes that quark mass effects 
are quite small for the radiatively corrected normalized angular decay distribution even if 
they are important for the radiatively corrected decay functions. 

When the polar angle 9 is measured w.r.t. the antiquark direction (z || q) (system II), 
the relevant expressions for the helicity structure functions can be obtained from those in 
system I by the exchange \i\ p,2, and H±± H TT and H 00 -h- H 00 . One has 

H ^ T (pi, /i 2 ) = H l ±± (p 2} fn), Hw(pi, P>2) = #oo(>2, fix)- (34) 

Because of the mass asymmetry of the polarized decay functions and because one is ex- 
changing m c "B- m&, the quark mass effects are more pronounced in system II. In fact, one 
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W + -> cb 


Z — > cc 


Z -> bb 




system I (z c) 


system II (z 6) 


sysiem i yz \\ cj 


sysiem i yz \\ u) 




—0.811 


—0.811 


+0.750 


+0.746 


c f 




-0.741 


-0.741 


+0.697 


+0.697 




-0.743 


-0.752 


_i_0 7D1 






+0.228 


-0.228 


+0.076 


+0.106 


A F b 


+0.224 


-0.224 


+0.073 


+0.102 




+0.225 


-0.226 


+0.073 


+0.103 




+0.281 


-0.281 


-0.101 


-0.142 


cos 9 


cxtr 


+0.304 


-0.304 


-0.105 


-0.147 




+0.302 


-0.301 


-0.105 


-0.146 



Table 1: The parameters Cf, Aps and cos 6? characterizing the normalized polar angle 

extr 

decay distribution of the cascade decays t — > b + lT r+ (t)(— > cb) and and the produc- 
tion/decay process e + e _ — > Z(t) — >■ 66, cc. We use a three-tiered notation where the top 
entries are for (Born; fa ^ 0), the middle entry is for (0(a s ); fa = 0), and the bottom 
entry is for (0(a s ); fa ^ 0). 

now has 

#?+ = 8N c q 2 [0.996+ ... + ^(1 + 1.806 + ...)], 



67T 



H, 



00 



II 



8N c q 2 [0.002 
8iV c g 2 [0+ .. 



-^(4-1.051+ ...)], 

67T 



+ — (1 - 0.284 + 

67T 



) 



(35) 



The NLO quark mass corrections in system II can be seen to be approximately two-and-a 
half times larger than those in system I. The largest mass correction now occurs for H+ + . 
Numerically one obtains 

_ q f 0.001 1 q 

W l \6) = -(1 + cosfl) 2 0.019 +-(l-cos 
8 y 0.016 J 8 

Again one observes that the large quark mass effects seen in the polarized decay functions 
do not feed through to the angular decay distribution. Of course, the reason is the same 



f 0.305 ) 




0.694 ) 




\ 0.318 \ 


+ |jsin 2 £<j 


0.663 \. 


(36) 


L 0.317 J 




0.667 J 
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(see Table 1). 

max 



as explained after Eq. (1291) . Except for the slightly enhanced NLO quark mass effects, the 
distribution ( 136]) is just a reflection of Eq. (|30|) at the line cos 9 = 0, i.e. A l FB ~ —A l p B . 

Larger quark mass effects can also be seen in the convexity parameter cj, in the forward- 
backward asymmetry A l p B and in the position of the maximum cos 9 
The larger size of the NLO finite mass corrections to the global measures can be inferred 
by listing approximate formulas similar to those in Eqs. ( 13TH331) . One has 

c l j = 7(1 -3a») (1-^(12 -3.153)) = -0.771, (37) 

3 4 V &7T ' 

A FB = - 7 (p ++ -p__)(l-^(6- 1.619)) = 0.222, (38) 
4 K on ' 

u P++-P— / 1 + ^£ (6 _ L534) \ = o.287. (39) 



cos 9 

1 - 3poo v 6tt 

One can see that the NLO finite mass effects make up ~ 25% of the leading NLO terms as 
compared to the ~ 10% in system I (see Eqs. ( I3TH331) ) . Again the approximate numbers 
listed in Eqs. (I37LI391) can be seen to only deviate by small amounts from the exact numbers 
listed in Table 1 (column 2; tier 3). Again one concludes that, even though quark mass 
effects are larger in system II, the bulk of the radiative corrections still come from the NLO 
leading terms cx 12a s /(67r) and oc 6a s /(67r). 

Up to this point we have only considered the decay W + — > q\q2- The charge conjugated 
decay W~ — > q\q2 is related to W + — > q\(j2 by CP-invariance. The corresponding helicity 
structure functions H mm (W~ — > q^) can be obtained via the relations 

H mm (W~ q x q 2 \ Hi, H2] z \\ q 2 ) = H mm (W + ->• q x q 2 \ /i 2 , P\\ z || q{), (40) 

where the z axis for the decay W~ — > q\q2 lies along the quark direction (z \\ q 2 )- 
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6 Polarized Z decays into massive quark pairs 

In the SM the Born term matrix element for the decay Z(m) — > qq with spin quantum 
numbers m is given by (see e.g. Ref. [5]) 



9z - - 
M z {m) = -i—u(q) [Vfj,, - a f ^^)v(q)e^(m) 



where g% = 8G F M§/V2, and where, in the SM, one has 

cif = 1 for u, c, t, 



v f 



1 - - sin 2 & w , 
4 

-1 + - sin 2 B w , 



for d, s, b. 



Similar to Eq. (fTTl) we define reduced amplitudes by writing 

Mzijn) = u(q) (v f 7 M - 0/7^75) v(q)^{m). 
As in Eq. ( fl"8l . the LO polarized decay functions H mm are calculated according to 

H mm = N c J2 M z {m)M z {m). 

quark spins 



(41) 

(42) 
(43) 

(44) 

(45) 



Compared to the charged current case the relative weights of the vector (V) and axial 
vector current (A) contributions are no longer simple and it is more convenient to switch 
to a notation in terms of the VV, AA and V A = AV contributions. Again we make use of 
the analytical results in Ref. [3] (or those in Ref. [13]) which we expand up to O(fi). One 
obtains 

_2 



VV/AA 
U 



H 



VA/AV 



H 



VV/AA 



±N c q 2 



4iV c g 2 



l-2/t±2/t ... + ^(2 + 2ttVm+(6 

D7T V 



57T 



)// 



-12/i In fj, + 2/i In 2 \i ± 24/a (1 + In /t) 



1 - 2/i . . . + ^ (32^ + (56 + - 8/1 In /1 



+2/iln /t . . . 



/i±/x... + ^(4-27r 2 v ^-(3^ 
6n v 3 

,2 



7T 



)/' 



-24/iln/z-2/iln>±6/i(5 + 21n/i) . 



(46) 
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Note that, in the zero mass limit, there are no a s corrections to the parity- violating struc- 



ture function Hp A , as noted before in Ref. [HI 115] - 



As in system I of the charged current case we evaluate the polarized decay functions 
H±± and H 00 in a system where the quark lies along the z direction [z \\ q). One has 



H i± = \(v)H^ + a)H^ T 2v f a f HV) 



ffJo = v)Hr + a)Ht A . (47) 

As a check on Eq. ( 14"T|) one can set Vf = a/ = 1 ort)/ = a/ = —1 and one will then recover 
the charged current results Eq. (fT9|) (system I) with ^1 — ^2 — When the z axis is taken 
to be along the antiquark direction (system II; z \\ q), there will be no change for H 00 but 
one needs to exchange H±± <->■ H TT (or, Hp A <->■ -Hp A ) in Eq. (j4"T|) . 
For the sum of the three polarized decay functions one obtains 



#[/+L = Hu+l = H++ + H oo + H -- = v ) h u+l + 0/ 
= 8iV c g 2 {^[J + // ... + ^(3 + 42// + ... 

I 1 L 2 D7T 



-2/i . .. + ^ (3 - 12//-36// In /i + ...)]}. (48) 

07T J J 



-2 

In order to make contact with the total decay rate T(Z — > gg), we define polarized decay 
rates by 

r = G F \p\ 

1D7TV2 

where \p\ = — 4/x. For the unpolarized decay rate one then obtains 

nz gg") = l(r ++ + r 00 + r__) = J^^W (so) 

Returning to Eq. ( I46p . we write down our numerical results for the decay Z — > bb. In 
system I one has 

— Cv 

H^ v (bb) = W c q 2 [l +-1(2 + 1.363 + ...)], 
H AA (bb) = 4N c q 2 [ 0.989 + ^(2 + 2.013 + •••)], 
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Hl A (bb) = 4:N c q 2 [0.994 + 77— (0 + 2.199 



HY V '(66) = 4iV c g 2 [0.006 + ^(4- 1.131 + ... 
i^(66) = 4iV c g 2 [0 (4 - 0.905 + ...)], (51) 

07T 



or, using Eq. (j4"7]l 



(66) = 4iV c g 2 [ 0.046 + ^(1.479 - 0.189 + ...)], 

67T 

#*(66) = 4iV c g 2 [ 0.003 + ^(5.916 - 1.447 +...)], 

— 01 
H l __(bb) = AN c q 2 \ 1.422 + -±(1.479 + 2.855 + ...)]• (52) 

6n 

The NLO quark mass effects in the polarized decay functions can be seen to be large. The 
largest NLO quark mass correction arises in the polarized decay function H l __(bb) where 
the mass correction amounts to a O(200%) effect compared to the leading NLO term. 
For the unpolarized decay function Hjj+l one obtains 

H u+L (bb) = AN c q 2 [ 1.471 + -± (8.874 + 1.219 + ...)] (53) 

67T 

The O(nb) quark mass and radiative corrections effects increase the zero mass Born term 
decay rate by 3.7%, i.e. one has (we take a s {m 2 z ) = 0.115) 

T( NLO; Z -> 66, 0(ji b )) = 1.037 • T( Born; Z -> 66, fi b = 0) . (54) 

Next we write down the angular decay distribution for Z(f) — > qq with Z polarization 
obtained from the production process e + e _ — > Z{\). In terms of the VV, AA and VA 
structure functions in Eq. ( 1461) one has 

W{9) = jj(l + cos 2 £) (v 2 f H™ + a 2 H^) 

3 2v e a 2v f a f 

4 v 2 + a 2 

+ 3 -sin 2 6 (v)Hr + a)H£ A ). (55) 

Note that the electroweak parameters v e and a e do not appear in the first and last row of 
Eq. ( 155]) because we are using normalized density matrix elements such that p ++ + p = 1. 
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Figure 2: Normalized angular decay distribution W(9) = W(9)/W for Z — > bb at LO (full 
line) and NLO (dotted line), including 0(m 2 ) finite quark mass contributions. 



The normalized decay distribution is obtained from Eq. fl55l) through W(9) 
whereW=(v 2 f HVV L + a 2 f H^ L ). 

Numerically one obtains [mz = 91.188 GeV, a s (m 2 z ) = 0.115) 



W{9)/W 



W\9) (bb) = -(1 + cosi 



0.570 
0.556 
0.559 



+ -(l-cos#) 5 



f 0.428 




0.002 ) 


\ 0.420 


j + 5sin^| 


0.024 \ 


l> 0.421 




0.020 > 



Eq. (|56l) describes an upward bent parabola with unit area and a minimum slightly 
displaced to the left of cos 9 = (see Fig. Final state radiative corrections have a 
0((1 — 2)%) effect (on an absolute scale) on the coefficient functions of the angular decay 
distribution while NLO quark mass effects contribute only at the per mill level. In Fig. [2]we 
compare the O(fi) angular decay distributions for the Born term case and the 0(a s ) case. 
The NLO corrections can be seen to make the angular decay distribution flatter. In fact, 
the convexity parameter is reduced by 5.36% through the radiative corrections as Table 1 
shows. The decay distribution is weighted towards the forward hemisphere such that A FB 
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is positive. It is barely visible that the radiative corrections reduce the forward-backward 
asymmetry. This is born out in Table 1 where one finds a 2.73% reduction in Afb- The 
minimum is slightly shifted to the left. Quantitatively, this amounts to a 2.78% effect as 
Table 1 shows. 

A rough description of the global effects of the radiative corrections on the decay distri- 
bution can again be obtained by approximate formulas using the same set of approximations 
as in Eq. ( l3~l~H33l) . One now obtains 

3, 



cos 6 



c\bb) 
A l FB (bb) 

L(») 



3p 



1 



QlT 



12 -2.93) 



3 

(P+ 



/ \ 2v f a f 



i 



= 0.708, (57) 
— (6- 1.37)) = 0.103, (58) 

07T / 



a 



P- 



-0.145. 



(59) 



1 - 3poo 

The numerical values obtained from the approximate formulas are quite close to the relevant 



numbers in Table 1 (column 4; tier 3) indicating that the approximation is quite good. The 
NLO finite mass effects in Eqs. f l57H59|) can be seen to reduce the respective leading NLO 
term by ~ 25%. This is in accordance with the numbers in Table 1 and similar to what 
happens in the decay W + — » cb (system II). 

The case Z — > cc has to be treated separately since, apart from the quark mass effects, 
one now has to use the electroweak coupling coefficients appropriate for up-type quarks. 
The numerical results are 



H V V V (cc) 



H* A (cc) 



a. 



VA 



H 



H V L V (cc) 



H£ A (cc) 



AN c q 2 [l +-^(2 + 0.367 + 
4N c q 2 [0.999 + — (2 + 0.460 + 

07T 

AN c q 2 [0.999 + ^ (0 + 0.599 + 



s 



4N c q 2 [ 0.001 + — (4 - 0.344 + 
6ir 



4N c q 2 [0 +-^(4 - 0.307 + 

D7T 



(60) 
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or, using Eq. (147"]) 



HlJcc) = 4N c q 2 [0.189 + -1(1.148 - 0.027 + ...)], 
FL(cc) = AN c q 2 [ 0.000 + ^(4.590 - 0.358 + ...)], 
Hl_(cc) = AN c q 2 \ 0.957 +^(1.148 + 0.488 + ...)]• (61) 

D7T 

For the unpolarized decay function Hu+l one obtains 

H u+L {cc) = AN c q 2 [ 1.146 + (6.886 + 0.103 + ...)] (62) 

on 

As expected, the mass corrections in the (cc) case are smaller than those in the (66) case. 
According to (|62|) the 0(fi c ) quark mass and the radiative corrections effects increase the 
zero mass Born term decay rate by 3.7%, i.e. one has 

T( NLO; Z cc, 0(> c )) = 1.037 ■ T( Born; Z cc, /x c = 0) . (63) 

For the normalized angular decay distribution one finds 

_ 3 f 0.551 1 q 

^(^(cc) = -(l + cos^) 2 0.537 +-(l-cos£) 2 
8 I 0.538 J 8 

The coefficients of the distribution are quite similar to those of the (66) case Eq. f )56|) 
only that the NLO finite mass effects are smaller than those in the (66) case. Apart from the 
small finite mass effects the shape of the angular decay distribution is also affected by the 
difference of the factor 2v fdjj (v 2 +a 2 ) in the (66) and (cc) cases which are given by 0.94 and 
0.67, respectively. This is relevant for the values of Afb and c/ (see Eqs. ( 158]) and ( 159]) ). As 
the relevant numbers in Table 1 show the Z — )■ cc forward-backward asymmetry is reduced 
by 28.9% relative to the (66) case. Also the minimum of the cc distribution is moved to 
the right by 28.4% going from the (66) to the (cc) case (see Table 1). 



r 0.449 




• 0.000 i 


<^ 0.439 




0.024 \ 


l> 0.440 • 




■ 0.022 ^ 
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7 Summary and Conclusions 

We have presented 0(a s ) results for the polarized decay functions that describe the decay of 
polarized (W, Z) bosons into massive quark-ant iquark pairs. NLO quark mass corrections 
to the polarized decay functions have been found to be quite large. They can be as large 
as 200% of the leading NLO mass term. However, these large NLO quark mass effects 
do not feed through to the normalized angular decay distributions. The large NLO quark 
mass effects disappear when one is dividing out the total decay function in the normalized 
angular decay distribution. 

We have combined these results with information on the spin density elements of the 
(W, Z) bosons in the two sample production processes t — > b + W + and e + e _ — > Z to write 
down explicit analytical and numerical forms of the polar angle decay distributions. The 
0(a s ) corrections to the polarized decay functions result in 0((1 — 3)%) absolute changes 
in the angular coefficients of the polar angle decay distributions where the bulk of the 
NLO corrections come from the leading mass term. Quark mass corrections are small but, 
depending on the required accuracy, are non- negligible. 

The radiative corrections make the angular decay distributions flatter for both W and 
Z decays where the main effect comes from the leading mass term in the final state NLO 
contribution. As a measure of the flatness of the decay distribution we have used the 
convexity parameter given by the second derivative of the decay distribution. For zero 
mass quarks the final state radiative corrections reduce the convexity parameters of the 
angular decay distributions by ~ (6 — 8)%. Depending on the particular case under study, 
NLO quark mass effects reduce this value to ~ (5 — 7)%. The forward-backward asymmetry 
is reduced by ~ 1.5% and ~ 3.5% for W and Z decays, respectively, where NLO finite 
mass effects have reduced these shifts by a small amount. The radiative corrections shift 
the position of the maximum or minimum of the decay distributions away from zero by 
~ 7% and ~ 3% for W and Z decays, respectively, where again NLO finite mass effects 
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reduce these shifts by small amounts. The total width of the decay W + — > cb is increased 
by 3.8% by NLO effects where NLO nonzero quark mass effects account for 7.8% of the 
increase. The total widths of Z — > bb and Z — > cc are increased by 3.7%, where NLO mass 
effects account for 13% and 1.6% of the increase, respectively. 

In this paper we have only discussed polar decay correlations which probe the diag- 
onal spin density matrix elements of the production processes. If one wants to probe in 
addition the nondiagonal spin density matrix of the gauge bosons, one needs to involve 
in addition azimuthal correlations between the momenta of the particles involved in the 
production/decay process. Upon azimuthal averaging one would recover the results of the 
present paper (see e.g. [10J). A comprehensive discussion of the azimuthal correlations 
would form the subject of a separate publication. 
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